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ABSTRACT 

The plateaus observed in about one half of the early X-ray afterglows are the most puzzling feature in 
gamma-ray bursts (GRBs) detected by Swift. By analyzing the temporal and spectral indices of a large 
X-ray plateau sample, we find that 55% can be explained by external, forward shock synchrotron emis- 
sion produced by a relativistic ejecta coasting in a p oc r~ 2 , wind-like medium; no energy injection into 
the shock is needed. After the ejecta collects enough medium and transitions to the adiabatic, deceler- 
ating blastwave phase, it produces the post-plateau decay. For those bursts consistent with this model, 
we find an upper-limit for the initial Lorentz factor of the ejecta, Tq < 46 (e e /0.1) _0 ' 24 (es/0.01) ' 17 ; 



the isotropic equivalent total ejecta energy is Ei s 



10 53 (e e /0.1)- 1 - 3 (e s /0.01)- - 09 (t 6 /10 4 s) erg 



where e e and eb are the fractions of the total energy at the shock downstream that are carried by 
electrons and the magnetic field, respectively, and tb is the end of the plateau. Our finding supports 
Wolf-Rayet stars as the progenitor stars of some GRBs. It raises intriguing questions about the origin 
of an intermediatc-Fo ejecta, which we speculate is connected to the GRB jet emergence from its host 
star. For the remaining 45% of the sample, the post-plateau decline is too rapid to be explained in 
the coasting-in-wind model, and energy injection appears to be required. 

Subject headings: gamma-ray burst: general - radiation mechanisms: non-thermal - relativistic pro- 
cesses - shock waves - X-rays: bursts 



1. INTRODUCTION 

Among the new features in early X-ray afterglows re- 
cently discovered by Swift, the X-ray plateau is the most 
troublesome to explain; see Zhang 2007 for a review. 
There is a long list of proposed explanations, including a 
prior emission model (Ioka et al. 2006; Yamazaki 2009), 
an evolving microphysical parameter model (Panaitescu 
et al. 2006; Ioka et al. 2006), a two-component (Granot 
et al. 2006) or multi-component (Toma et al. 2006) jet 
model, a slow energy transfer from the magnetized ejecta 
to the ambient medium (Kobayashi & Zhang 2007), a 
reverse-shock dominated afterglow model (Uhm & Be- 
loborodov 2007; Genet et al. 2007), a circumburst-dust- 
scattcrcd echo of prompt X-rays (Shao & Dai 2007; but 
see Shen et al. 2009), and a scattered internal shock 
(Shen, Barniol Duran & Kumar 2008) or external shock 
emission (Panaitescu 2008) model. Most fail to explain 
all the properties of the X-ray plateau. The leading 
model invokes energy injection to refresh the external 
shock. In this scenario the injected energy comes either 
from a long- lasting central engine activity (e.g., Dai & Lu 
1998a, b; Zhang & Meszaros 2001; Dai 2004; Yu & Dai 
2007) or from a slower portion of an outflow with a broad 
initial Lorentz factor (LF) distribution (e.g., Granot & 
Kumar 2006). However, the standard refreshed shock 
model faces several serious issues. First, the flat slope 
(~ t~ - 5 ) and late ending time (~ 10 4 s) of the plateau 
feature imply that the total energy carried by the later or 
the slower ejecta is much larger than that of the prompt, 
fast component which gives rise to 7-rays. Second, the 
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abrupt end to the plateau phase poses a serious theo- 
retical challenge to models for the origin of the injected 
energy. The late-activity scenario also suffers from the 
fact that prompt 7-rays, which are comparable in energy 
to the late-time injection, emerge so much earlier. The 
implied radiative efficiency can exceed 90% (e.g., Ioka et 
al. 2006; Zhang et al. 2007), much higher than the ~ 1% 
efficiencies attained by the popular internal-shock model. 
We note that if the energy injection is made by an pair 
dominated outflow, such as a dissipated wind from a mil- 
lisecond magnetar (Dai 2004; Yu & Dai 2007; Yu, Liu & 
Dai 2007; Mao et al. 2010), the wind-decelerating rela- 
tivistic reverse shock radiates more efficiently than in the 
case of a baryon-dominated injecting outflow and might 
even dominate over the forward shock emission; this al- 
leviates the energy budget and the required efficiency. In 
that model, the X-ray plateau ending time corresponds 
to the spin-down time of the magnetar. 

Some afterglow features under the fireball shock model 
in a p cx r~ k medium, where k is a constant, have been 
described earlier by Dai & Lu (1998c) and Meszaros, Rees 
& Wijers (1998). Chevalier & Li (1999, 2000) have stud- 
ied specific behavior of afterglow due to the blastwave 
expanding into the wind of a Wolf-Rayet star, corre- 
sponding to k = 2. Panaitescu & Kumar (2000) stud- 
ied the afterglow light curves both in a constant den- 
sity {k = 0) interstellar medium (ISM) and in a wind 
medium. Models of late-time (t > 10 4 s) afterglow data 
(e.g., Panaitescu & Kumar 2002; Schulze et al. 2011; 
Oates et al. 2011) have inferred that afterglows consis- 
tent with ISM are more numerous than those consistent 
with a wind medium. 

In the external shock model for gamma-ray burst 
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(GRB) afterglows, relativistic ejecta coasts freely before 
it collects 1 /Tj times its rest mass from the circum- burst 
medium and decelerates, where Tj is its initial LF. In 
the past, afterglow studies have concentrated mainly on 
the decelerating phase, in which the flatness of the X-ray 
plateau decay slope is difficult to explain. In this paper 
we look at the early, coasting phase. If the medium den- 
sity is constant with radius, the light curve in this phase 
should rise as t 2 ~ 3 , simply because the total number of 
emitting particles increases steeply with time. A feature 
consistent with such a rise was only detected for a few 
cases (e.g., Molinari et al. 2007). Such a rise may of- 
ten be difficult to detect if it finishes so early that it is 
buried below the prompt emission tail, which is dimming 
but still brighter than the rising afterglow. 

In a wind medium, however, the external shock light 
curve in the coasting phase is flat or slowly decaying, if 
the observing frequency is above both the synchrotron 
injection frequency and the cooling frequency (see $21 
and Table [1]). Waxman (2004) discussed a mildly rel- 
ativistic shock coasting in the wind for explaining the X- 
ray plateau in GRB 980425 / SN I998bw (Kouveliotou 
et al. 2004) long before the Swift era. In this paper 
we investigate systematically the scenario in which the 
SW/i-observed X-ray plateau is produced in the coasting 
phase of the ejecta plunging into a wind medium, and 
the 'normal' decay following the plateau corresponds to 
the subsequent decelerating phase of the ejecta in the 
same wind medium; we refer to this as the "coasting-in- 
wind" model. We assume the synchrotron emission of 
the external forward shock is the dominating emission 
component. 

We start with reviewing in Sj2]the dynamics and radia- 
tion properties of the external forward shock in a general 
medium, for both the coasting and decelerating phases of 
the ejecta. Then in $3] we model a large sample of X-ray 
plateau data with the coasting-in-wind model, and in J2] 
we derive constraints on the initial LF and the isotropic 
equivalent energy of the ejecta. Summary of the results 
and discussion of the implications are given in <j5j 

2. EMISSION FROM EXTERNAL FORWARD SHOCK IN A 
GENERALIZED MEDIUM 

In this work we use the simplest self-consistent assump- 
tions and approximations, including assumptions that all 
electrons are accelerated into a single power-law distri- 
bution, and that the magnetic and electron energy den- 
sities are fixed fractions of the post-shock energy density. 
Simulations of relativistic shock acceleration (Spitkovsky 
2008; Sironi & Spitkovsky 2009, 2011) show that only a 
small fraction (~ 1%) of downstream electrons are ac- 
celerated in a power-law distribution carrying ~ 10% of 
downstream energy. This can be modeled using a fixed 
e e (defined below), unless in reality either the number 
fraction or the energy fraction of non-thermal electrons 
changes strongly with the shock LF (e.g., see £j3.2[) . The 
standard afterglow modeling also assumes that electron 
Larmor radii are smaller than the coherence length of 
the magnetic field so that synchrotron radiation applies; 
otherwise the emission is in the so-called 'jitter' regime 
(Medvedev 2000). 

After crossing the shock, electrons are accelerated into 



a power-law energy distribution with an index p. Those 
non-thermal electrons carry a small fraction, e e , of the 
internal energy downstream. The minimum electron en- 
ergy in this power law is 

!f p e e T(j, e m p /m e , for 2 < p, 

(f p e e r-flf 2 fi e m p /m^J , for I < p < 2, 

where p e is the number of nucleons per electron, f p — 
\p— 2\/(p— 1) derives from integrating the number density 
over the power law (e.g., Sari et al. 1998), and for the 
second case 7m = \/3e/ (ctB) is the maximum electron 
energy in the power law (Dai & Cheng 2001). Here <tt is 
the Thomson cross section, e is the electron charge and B 

is the magnetic field, whose strength B = y / 32neBT 2 pc 2 
is determined by the fraction, es, of the downstream en- 
ergy which goes to the magnetic field. The accelerated 
electrons emit synchrotron radiation. Another critical 
electron energy is the cooling energy above which elec- 
trons lose significant energy in the dynamical time: 

_ 67rm e c(l + z) 
lc ~ (1 + Y)a T TBH> [ ) 

where z is the GRB host red shift and Y is the syn- 
chrotron self Compton (SSC) parameter; we neglect the 
SSC loss, i.e., Y = 0. Above 7 C , cooling modifies the 
power law slope of electrons distribution. Radiative 
losses are necessarily important when p < 2, because in 
this case the electron energy is concentrated at very high 
energies above 7 C ; conversely when p > 2 synchrotron 
losses are important only when j c < j m . 

For the dynamical model, we consider a relativistic 
ejecta of isotropic equivalent energy Ei SO and LF Tj 
plunging into the circum-burst medium. The medium 
density profile is generally assumed to be a power law 
p = Ar~ k , where for the steady stellar wind case, 
A k=2 = M w /(4nv w ) = 5 x 10 n A* g cm' 1 (Chevalier 
& Li 1999); the fiducial value A* = 1 would arise for a 
wind of mass loss rate M w = 10 _5 M Q yr _1 and speed 
v w = 10 3 km s _1 . In the case of a uniform density ('ISM') 
medium, A^-a = p. 

A pair of shocks is produced, with the forward shock 
(LF F s h) moving into the medium and the reverse shock 
into the ejecta. The shocked fluid region comprises 
the forward-shocked medium and reverse-shocked ejecta, 
which have approximately the same LF T(t) w r 8 h/v2 
but are separated by a contact discontinuity. T(t) is 
given by pressure balance at the discontinuity as T — 
Tj/(1 + 2r i /a 1 / 2 ) 1 / 2 , where a = p'j/p{r) is the rest- 
frame density ratio between the ejecta and the ambient 
medium (e.g., Panaitescu & Kumar 2004b). In the limit 
of a» F 2 j, T ~ Tj. It can be shown that, before the de- 
celeration, i.e., the reverse shock passage, of the ejecta, 
the dense shell condition a ^> T 2 is satisfiedQ In the 
coasting-in-wind model Tj takes a single value rather 

1 The ambient density is p(r) = (3 — k)M(r)/(4nr 3 ), and the 
rest- frame jet density is p'^ = Mj / (47rr 2 A^. ) , where Mj is the 
isotropic ejecta mass and A^. is the ejecta width in the rest frame. 
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than a distribution, so that T(t) has a definite initial 
value To ~ Tj. 

Right after the reverse shock passage, about half of the 
initial ejecta total energy is transferred to the shocked 
medium and ejecta. The combined shocked region enters 
the Blandford & McKee (1976) self-similar solution, in 
which T(t) declines. For convenience, we denote phase 1 
as the period in which the reverse shock is crossing the 
ejecta, and phase 2 the later period of blastwave deceler- 
ation. In phase 1, since usually the reverse shock is too 
weak to account for the X-ray emission, we consider the 
forward-shocked region only whenever radiation proper- 
ties are concerned. 

In phase 2 the shock compresses the density and en- 
hances the energy per unit mass both by a factor of ~ T; 
energy conservation requires (Cohen et al. 1998; also see 
Wu et al. 2005 for similar treatment of radiative loss) 



- sh' 



17-4fc + e 
where, to account for synchrotron losses, 



0, 



if p > 2 and 7„ 
otherwise 



< 7c 



(3) 



(4) 



is the fraction of post-shock energy, if any, radiated away; 
approximation (j4|) matches the more complicated exact 
solution by Cohen et al (1998) to within 5%. If p > 2 
and 7 m < j c so that E- lso is conserved in phase 2, T cx 
r -(3-k)/2 K t -(3-fc)/(8-2k). oth e rw i S e T 2 cx r 

(Cohen et al. 1998) 

(l + e) 2 +3(l + e)(4-fc)-4 

m = 

3-e 

which in the adiabatic limit becomes m e= o — 3 — k. In 
any phase of continuous motion with T 2 oc r~ m , the ob- 
server's time t is given by 

1 + z r 1 + z r 

(6) 



where 



(5) 



t 



1 + m ctT 2 c 



4-kc t r 2 c 



where the arrow evaluates m — » (3 — k) for a deceler- 
ating, adiabatic blastwave. For adiabatic blastwaves, 
we find r cx r~ x / 2 oc i -1 / 4 in a wind medium where 
t = (1 + z)r/(2c t r 2 c), and T oc r~ 3 / 2 oc t~ 3 / 8 in a uni- 
form medium where t = (l + ^)r/(4ctr 2 c). The pref actor 
ct = 4 if most of the observed emission is from a thin layer 
of shocked fluid right behind the shock front and on the 
line-of-sight axis (Sari 1997; Dai & Lu 1998c; Chevalier 
& Li 2000). Emission from off axis and from the fluid 
further downstream arrives somewhat later, so the effec- 
tive Ct is somewhat smaller (Waxman 1997; Sari 1998; 
Panaitescu & Mesazaros 1998). We adopt c t = 2. 

The observed synchrotron frequency for electrons with 
energy 7, averaged over an isotropic distribution of 
the electron's pitch angle, is v = 3r7 2 e-B/[16m e c(l + 

For a late deceleration time, which is appropriate for our coasting- 
in-wind model, the ejecta width is determined by radial spreading, 
so Aj. ~ r/Tj. Thus a = p^/p(r) ~ TjMj/M(r). Prior to decel- 
eration, Mj/Fj S> M(r) (cf. Equation [3]l , so o > T?; however 
a ~ at deceleration. 



z)\. The critical frequencies corresponding to 7™, 
and 7 C are v m and v Cl respectively. The peak syn- 
chrotron specific power for a single electron is P^max = 
8m e c 2 aT^B /(9we), independent of the electron energy. 
The peak flux density in the observed F v spectrum 
is i^max = A^maxU + z)/{4^eP L ), where N = 
4:TrAr 3 ~ k /[(3 — k)m p ] is the total number of electrons in 
the swept-up medium, and is the luminosity distance. 

Given the above formulae, we can find the following 
time evolution laws for the case of p > 2: 

u m oct- k / 2 7 ^oci 3fe / 2 - 2 , i7,, max cxi 3 - 3fc / 2 (7) 

for phase 1, and 



v m oc t 2 <™+ 1 ) , ^ c oct 2 (™+!), F v 



oc t 



(8) 



for phase 2. For 1 < p < 2, the only difference is on 

-k _ m(p + 2) + k 

v m : v m oc 1 2 (p-!) in phase 1, and v m oc t 2 ( m +i)( P -i) m 
phase 2. Note that for wind medium, v m decreases and 
v c increases monotonically with time in both phases. 

Following from the above evolution laws the observed 
flux density F v evolves according to F v (t) oc v~Pt~ a 
(e.g., Sari, Piran & Narayan 1998). The spectral slope 
/3 is determined by p and the relative locations of v 1 v m 
and v c . Therefore, a specific numerical relation between 
a and (3, the so-called "closure relation" , exists depend- 
ing on which type of medium and which dynamical phase 
the external shock is in. We refer the reader to Zhang 
& Meszaros (2004) for a summary of these relations in 
the decelerating phase for both uniform and wind media. 
For demonstration purpose, we list in Table[T]the closure 
relations for both the coasting phase and the decelerating 
phase in a k — 2 wind medium. 

3. SAMPLE AND DATA ANALYSIS 

We use the data from a sample of 53 X-ray plateaus 
compiled by Liang et al. (2007). This sample is provided 
with fitted spectral and temporal slopes of the plateau 
component, /3j and a±, and of the normal decay compo- 
nent following the plateau, /?2 and 0:2, and the time of the 
transition between the two components tt>. We omit four 
bursts whose phase 2 decay slopes are too steep ~ 3 
for GRB 060413, 060522 and 060607A; a 2 « 9 for GRB 
070110) to be explained by external shocks; for them the 
best explanation is that some internal dissipation of a 
delayed central engine ejecta is responsible for both the 
abrupt step-like decline in phase 2 and the prior plateau 
(e.g., Troja et al. 2007; Lyons et al. 2010). This leaves 
us a sample of 49 bursts. 

Figure [T] compares the observed values of a.\ and (3i 
against the predicted closure relations for different dy- 
namical, medium and spectral models. It can be seen 
that the coasting-in-wind model with max(j/ m , u c ) < vx 
most successfully accommodates the majority of data. 

Within this model the post-plateau decay corresponds 
to the decelerating phase (phase 2) of the blast wave in 
a wind medium; the break cannot be due to a change 
of spectral regime because there is no spectral evolution 
(#2 = Px) in all but one burst — GRB 061202 (e.g., 
Liang et al. 2007; Shen et al. 2009). The lack of spectral 
evolution also signals a requirement for any model, that 
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TABLE 1 

THE SPECTRAL AND TEMPORAL INDICES OF THE FORWARD SHOCK SYNCHROTRON EMISSION (Ft, OC V~@t~ a ) IN THE COASTING-IN-WIND MODEL 

FOR THE WIND DENSITY INDEX k = 2. 



Spectral regime 



a (p > 2) a (p < 2) 

phase 1 phase 2 phase 1 phase 2 



Slow cooling (v m < u c ) 
v < v m 

Um < U < U c 



_1 _1 n 1 5(2-p) 
3 3 u 3(p-l) 12(p+l) 
J-l p-1 3p-l 1 p+8 
_2 2 i 2 2 



Fast cooling (u c < Vm) 

v < u c 

Vc < V < V m 



111 1 2 

3 3j 3 3j j 



_2 2, .. 4 „ 2_ 



P+b 



max(i/ m , f c ) < 



~P p— 2 ;-ip — 2 

2 2 4_ 



Note. — Phases 1 and 2 are the coasting and decelerating phases of the blast wave, respectively. 



1 - 



- 




_1 . Coasting in Wind max(v m ,v c ) < v x 

— — - Coasting in Wind v m < v x < v c 
■ — ■ — ■ Decelerating in Wind v m < v x < v c 
Decelerating in ISM v < v (/ < v 

s m X c 

0.4 0.6 0.8 1 1.2 1.4 1.6 




Fig. 1. — X-ray decay index vs. spectral index during the plateau 
phase. The lines are the closure relations predicted by various 
models. The prediction from the coasting-in-wind model with 
max(y m ,^ c ) < v (thick solid line) is the most capable one that 
accounts the data. The predictions from the same model with 
slight variations in k (thin solid lines) are also plotted. The data 
in open circles (magenta color) are those whose post-plateau de- 
cays are inconsistent with the decelerating phase (phase 2) of the 
coasting-in-wind model (see Figure[2]and in the text). GRB 061202 
that shows a strong spectral evolution is plotted as a diamond. 

is, whichever spectral regime phase 1 is in, it remains in 
phase 2. We therefore expect a specific relation between 
«2 and /?2 depending on the relative locations of ux , v m 
and v c . 

Figure [2] plots the observed a 2 vs. fa with pre- 
dicted relations superimposed. For the decelerating 
phase in a wind medium, we see that the spectral regime 
max(^ m ,^ c ) < vx is the most consistent with the data, 
which also satisfies the requirement of zero spectral evo- 
lution. However, as in Figure [TJ the data in Figure [2] 
show broad scatter around the prediction of this model. 
While the scatter in the former figure could be attributed 
to a slight variation of the density index k from burst 
to burst, the scatter in the latter cannot — because oti 
has no ^-dependence in the max(!/ m ,v c ) < vx spectral 




Decelerating in Wind max(v m ,v c ) < v x - 

Decelerating in Wind v m <v x <v c 

— — ■ Decelerating in ISM v m <v x <v c 



0X,2 

Fig. 2. — X-ray temporal index vs. spectral index during the 
post-plateau phase. The solid lines are for the decelerating phase of 
the coasting-in-wind model. Note that for a decelerating blastwavc 
with max(f m , u c ) < v and p > 2, a is independent of k. The data 
in open circles (magenta color) are for the inconsistent sub-sample. 
The thick dotted line (blue color) is the prediction for phase 2 of the 
coasting-in-wind model considering the radiative loss of blastwave 
energy (requiring v c < v m when p > 2, and assuming e e = 0.3); it 
should be taken as an upper limit of ctx.2 for a given f3x,2 within 
this model because of the high e e value assumed. The exceptional 
burst that shows a spectral evolution is labeled. 

regime. This is a puzzle for all decelerating blast wave 
models that do not invoke late energy injection. 

In conclusion, we find that the coasting-in-wind model 
for the plateau, and its subsequent (adiabatic) decelerat- 
ing stage for the post-plateau decay, can account for 55% 
of the sample. We call these bursts the consistent sub- 
sample, and they are plotted as black data points in Fig- 
ures [T] and [2] GRB 061202 is consistent with the closure 
relation predictions of the coasting-in-wind model both 
during and after the plateau phase, but its strong soft- 
ening evolution from /?i = 1.15±0.07 to j3 2 = 2.55 ±0.44 
cannot be accounted for by a nominal spectral regime 
change v m < vx < v c — > max(v m , v c ) < fx, nor by any 
variant of the coasting-in-wind model; we therefore ex- 
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elude it from the consistent sub-sample. The remaining 
of the sample, i.e., those located above the coasting-in- 
wind model prediction by at least 2a in Figure [2] are 
marked as open circles (in magenta color) both in Fig- 
ure [2] and Figure [1] We call these the inconsistent sub- 
sample. We see in Figure Q] that the inconsistent sub- 
sample shows no significant difference from the consis- 
tent sub-sample in the a.\ vs. /3i scatter plot. Figure [3] 
shows the distribution of k derived from the consistent 
sub-sample in phase 1, which is centered around k = 2. 

At face value, this finding seems to contrast with pre- 
vious studies (e.g., Panaitescu & Kumar 2002; Schulze 
et al. 2011;Oates et al. 2011) of late-time (t > 10 4 s) 
afterglow data, which concluded that GRBs with k ~ 
predominate over those with k ~ 2. However, the dif- 
ference should not come as a surprise. Our consistent 
sub-sample includes only those bursts consistent with the 
max(f m , v c ) < vx spectral regime — as required by our 
modeling of the plateau data and by the lack of spectral 
evolution. Since k cannot be determined from the model 
prediction in phase 2 for this spectral regime (solid line 
in Figure [5]), we infer k instead from plateau data. Our 
inferred k values, based as they are on the assumption 
of a coasting ejecta, are model dependent: we would ob- 
tain very different values within a model involving late 
central engine activity. 

As for the inconsistent sub-sample, interpreting these 
under the v m < vx < v c spectral regime indicates that 
many are closer to k — than k — 2. Two bursts, 
050801 and 060714, appear in both the Schulze et al. 
(2011) sample and our consistent sub-sample; Schulze et 
al. classify both as k = 0. Such overlap is expected, 
as we describe a burst as 'consistent', and derive its k 
from the coasting-in-wind model, if it is within 2a of the 
max(^ m , v c ) < vx locus in phase 2. Bursts at the upper 
periphery of this sub-sample could also be interpreted 
assuming v m < vx < v c , and the derived k would be 
very different. Note also that our values for the two 
bursts' post-plateau indices, from Liang et al. (2007), 
differ slightly from those used by Schulze et al., and the 
latter authors treat a burst as consistent if it falls within 
3a of a model, whereas our threshold is 2a. 

Allowing for the radiative loss of blastwave energy in 
phase 2 can help to increase the size of the consistent 
sub-sample, but not by much. For instance, in the most 
radiative case — e e = 0.3, and v c < v m for p > 2 particu- 
larly — the predicted ax, 2 in Figure [2] is shifted upward 
by w 0.2. That freedom could at most increase the con- 
sistent sub-sample size to about 2/3 of the total. 

To account for the inconsistent sub-sample, one must 
invoke the energy refreshment in the decelerating blast 
wave. That is the topic of the next subsection. 

3.1. Refreshed shock models 

There are two primary scenarios in which the blast 
wave can be refr eshecQ. 

2 A third scenario was suggested by Kobayashi & Zhang (2007), 
in which the ejection is brief but is dominated by the Poynting 
flux. The magnetic energy in the ejecta is not transferred to the 
shocked ambient medium until the ceasing of the reverse shock. 
The delayed transfer of the magnetic energy serves as a varied 
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Fig. 3. — Distribution of the derived density power-law index k 
for the sub-sample that are consistent with the coasting-in-wind 
model in both its phase 1 and phase 2 — the consistent sub-sample. 
The top left horizontal bar shows the size of the typical error of k. 

1. The central engine activity is long lasting (Zhang 
& Meszaros 2001) with an outflow luminosity history 
L(t) oc t~ q . Since the LF of the outflow is always much 
larger than that of the blast wave, this is also the rate at 
which the total energy of the blast wave gets replenished, 
i.e., Ei so {t) oc t l ~ q . In order for the late energy injection 
to be important, q < 1. 

2. The central engine is briefly active, but after the 
prompt emission phase finishes, the ejecta has a wide 
range of LF Tj such that outer shells have larger Tj and 
moreover the ejecta mass (or even energy) is dominated 
by the inner, slower shells, i.e., the mass of ejecta is a 
power law in Tj (Rees & Meszaros 1998; Panaitescu, 
Meszaros & Rees 1998; Sari & Meszaros 2000; Ramirez- 
Ruiz, Merloni & Rees 2001): M(> Fj) cx TJ S . This 
extended ejecta distribution could be the outcome of a 
GRB jet breaking out of the progenitor star surface or the 
internal shocks during the prompt emission phase. The 
coasting-in-wind model can be considered the s — > oo 
limit of this class, in the sense that the ejecta have a 
single Tj. 

In both scenarios the blast wave energy is replenished 
at some certain rate, so as to cause a flattened plateau 
temporal slope, and tj, corresponds to the ending of en- 
ergy injection; the post-plateau decay is the standard 
decelerating blast wave phase. For any given ct\ and 
a medium type, there is a one-to-one correspondence 
relation between values of q and s for which the for- 
ward shock evolution is identical (although the reverse 
shock is not) , so it is not possible to distinguish between 
the two scenarios on the basis of X-ray data alone. We 
find that the inconsistent sub-sample can be accommo- 
dated within these models in either a wind or a uniform 

version of the two scenarios mentioned here for the refreshed shock 
model. However, the process of the magnetic energy transfer in 
the shocked region is poorly understood. Fan & Piran (2006) and 
Mimica et al. (2009) claim that the time scale for the transfer is at 
most several times the light crossing time of the ejecta, too short 
to account for the plateau feature. 
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Fig. 4. — Distribution of the derived q (top) and s (bottom) for 
the inconsistent sub-sample that require a refreshed shock for in- 
terpretation, where q is the late central engine luminosity temporal 
index and s is the catching-up ejecta mass distribution index (see 
section [XTJ . The spectral regime v m < ux < Vc is required by the 
data in phase 2. For a refreshed shock in a wind medium, the data 
require unreasonable conditions, i.e., q < and s < 0, for the two 
shock refreshment scenarios, respectively. Therefore, ISM is the 
only allowed medium type when a refreshed shock is considered. 



medium, so long as the observed band is below the cool- 
ing break (y m < v < v c ). The lack of spectral evolution 
implies that this is true in the plateau phase as well. 
With energy injection a new closure relation exists in 
phase 1 among ai, /3i and q (or s), within a given ambi- 
ent medium (Zhang et al. 2006). From this we can infer 
the required q (or s) value for each burst in the inconsis- 
tent sub-sample. The resultant distributions of q and s 
are shown in Figure |4j 

Figure[4]shows that within a wind medium the required 
parameters are unreasonable, in the sense that q < or 
s < for the two varieties of model. The inconsistent 
sub-sample can therefore be best accommodated within 
a refreshed-shock model (of either variety) operating in 
a uniform medium, and the cooling break must be above 
the observed band. Note, however, that in the post- 



plateau phase a small portion of these bursts appears to 
decelerate within a wind rather than a uniform medium 
(see Figure [2] the dashed line) ; further examination of 
these is warranted. The inferred ^-distribution peaks at 
q ~ 0.5 and the s-distribution peaks at s ~ 3. This 
means that in these two refreshed shock scenarios, re- 
spectively, the outflow total kinetic energy is dominated 
either by the late ejecta or by the slower, massive ejecta. 

3.2. Optical data 

A minority (13) of the original X-ray plateau sample in 
Liang et al. (2007) have simultaneous optical monitoring 
data; six of them are within our consistent sub-sample. 
Comparison of the optical light curve with the X-ray light 
curve shows diversity. Three (GRB 050801, 060714, and 
060729) show achromatic breaks around tb, while in the 
other three the breaks are chromatic. This roughly equal 
ratio of achromatic to chromatic bursts is about the same 
for the total 13 X-ray plateaus with optical data in Liang 
et al. (2007). Among the latter three, GRB 060210 
shows an optical break much earlier than tb (there is no 
optical coverage on % or later), GRB 060526 shows an 
optical break later than tb, and the optical light curve in 
GRB 050319 shows no break at all and is consistent with 
a single power law fit. Can our coasting-in-wind model 
explain the diverse optical behavior as well? 

The achromatic breaks seen in the former three bursts 
are consistent with the coasting-in-wind model, accord- 
ing to which the break at tb is of dynamical origin. In 
these three bursts the optical and X-ray temporal slopes 
are very similar, both before and after the break (see 
Figure 2 and Table 2 of Liang et al. 2007). This implies 
the same physical origin for emission in the two wave 
bands, which fits well to the coasting-in-wind model. 
This consistency also applies to the sub-sample compiled 
by Panaitescu & Vestrand (2011) whose coupled optical 
and X-ray light curves both show plateaus (their Figure 
4). 

An early optical break before tb, such as seen in GRB 
060210, can be interpreted as well in the coasting-in-wind 
model by the passage of v m or v c across the optical band 
during phase 1; for instance, the spectral regime change 
v < v m < v c — > v m < v Q < v c causes the slope steepen 
by Aa = (3p — l)/6 (see Table [IJ cf. observed Aa f=s 1 
in GRB 060210). Note that this change is consistent 
with the restriction that v m decreases and v c increases 
in the model. This pure spectral change scenario could 
also explain those bursts in which the optical shows a 
peak at the beginning of or during the X-ray plateau 
(e.g., GRB 060607A and 061121; see Liang et al. 2010 
and Panaitescu & Vestrand 2011). Similarly, a v m or v c 
passage over v Q in phase 2 can explain the optical break 
later than tb, such as seen in GRB 060526. 

The real difficulty lies in explaining the lack of an opti- 
cal break around % as seen in GRB 050319. If the end of 
X-ray plateau is of a dynamical origin, as postulated in 
the coasting-in-wind model, it is hard to avoid an achro- 
matic break at tb- Panaitescu et al. (2006) proposed 
that one solution is to let e e and es be functions of the 
blastwave LF T (also see Panaitescu & Kumar 2004a): 



oc T~ 



e B oc T~ 



(9) 
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We adopt this treatment here for the coasting-in-wind 
model, which means in phase 1 e e and e# are still con- 
stant, but they are free to evolve with time in phase 
2. Note that Panaitescu et al.'s scenario uses energy in- 
jection for the X-ray plateau while the coasting-in-wind 
model does not. In the following, we derive the condition 
on e and b for the lack of optical break at t b . 

We first calculate the X-ray decay slope in phase 2 in 
the presence of evolving e e and e^, since we know its 
spectral regime is max(v m ,y c ) < vx- 



1 b 

ax 2 = — I 

' 4 16 



(P ~ J ) (I - t - I 
2 2 8 



(10) 



Next we calculate the optical slope in phase 2. The ob- 
served slope a in those single power-law optical light 
curves usually lies in the range of 0.5 - 0.8. In phase 
1, this slope is best consistent with the spectral regime 
v m < v Q < v c for which a 0j i = (p — l)/2 (see Tabled]). 
This spectral regime should remain in phase 2 since there 
is no break in the light curve. Therefore, in the presence 
of evolving e e and , 



1 b 



(P-1) 



(11) 



Then, the equality a = a .i = a Q .2 poses the following 
condition for the lack of an optical break: 



a a 



(12) 



In the case of GRB 050319, a a R 
becomes 

36 + 4e = 12. 



0.5. The condition 



(13) 



Together with the observed a x .i ~ 1 for GRB 050319 
and Eq. (fTUfl . it gives e ps and b ps 4 as the condition 
for a single power-law decay in optical contemporaneous 
to a plateau-to-normal decay transition in X-rays. 

A few other bursts show the similar chromatic behav- 
ior as in GRB 050319 (e.g., Panaitescu et al. 2006; Oates 
et al. 2011). However, as pointed by Panaitescu et al. 
(2006), there is no obvious reason that the evolution of e e 
and cb has to satisfy such relation as in Eq. (fl3|) . This 
is the difficulty of this scenario. Alternatively, the chro- 
maticity in these bursts suggests that the X-ray and op- 
tical afterglow emission may arise from different origins, 
for instance, from a two-component outflow in which a 
low -Eiso, high Tj component decelerates very early and 
produces the optical emission with a single power law 
decay (the LF of this jet component has dropped signif- 
icantly from its initial Tj, which explains its dominance 
in the long- wavelength emission) , and a high E- lso , low Tj 
component decelerates later at tb and is responsible for 
the X-rays. 

4. CONSTRAINTS ON PARAMETERS OF THE 
COASTING-IN-WIND MODEL 

In the previous section we examined a large sample 
of X-ray plateaus and found the coasting-in-wind model 
with max(i/ m ,i/ c ) < vx to be consistent with the ma- 
jority of bursts. While this model is consistent with the 



data for all the sample in the plateau phase, it can ac- 
commodate only 55% of sample data in the post-plateau 
phase. In this section we aim to put constraints on the 
model parameters based on the consistent sub-sample. 

For this model to work, the deceleration time t dcc must 
equal t b \ t dcc is given by setting T(t dcc ) = T in Eq. 
([3]). Following from this, the blast wave isotropic en- 
ergy can be inferred. The synchrotron spectral regime 
max(i/ m , v c ) < vx must be justified for both phase 1 and 
2. Since with time v m decreases and v c increases mono- 
tonically during both phases (Eqs. [7] - [8]) , the spectral 
constraints are: v m (tx) < vx and v c (t2) < vx-, where 
t\ is the observed time of the earliest plateau data point 
and ti is that of the latest post-plateau decay data point. 
The last constraint is that the predicted flux density level 
has to match the observed one. Without losing general- 
ity, we choose to calculate the flux density at f&- Those 
constraints are summarized as follows. 

• ^dec = h ~ 10 4 s. 

• v m {t = t\ ~ 500 s) < vx ■ 

• v c (t = t 2 ~ 10 5 s) < v x - 

• F vx (i = t b ~ 10 4 s) - 1 /iJy. 

The formulae to calculate the model predictions are 



= 4 ' 5 X gl k x 10 52-10fc AT Z-2k 



17 -4k 



o 



ti, 



1 + z 



3-k 



erg. 

(14) 



fm(*l) = < 



'3 x 10 24 - 5fe x 6- fc / 2 /|ef4 /2r o~*^ 1/2 *r fc/2 

x (1 + zfl 2 - 1 Hz, for 2 < p, 
6.5 x 10 21 

x (4.8 x 10 2 ~ 5fe x 6- fc / 2 / 2 e 2 e ^ /2 A 1 /2)i/(p-i) 

x iy^^ 2 ^- 1 '^ + ^W^rr- 1 Hz, 
for 1 < p < 2, 

(15) 

v c (t 2 ) = 0.65 x [6(4- fc)] 3fe / 2 x 10 15fe - 10 e B 3/2 ^- 3 / 2 rf- 4 

, \ 3fe/2-2 
t2 \ 



1 + Z 



_ (3-fc)(3fc-4) 

(16) 



0.72xl0 3 x(3.1xl0 3 ) p ~ 2 x(6xl0 10 ) (P j 2) k 

(3-fc)(i+5)5 _1 £'i 8 

v fp-l £ p-l c V ,£+^ r (p+2)(2-fc/2) 
x Jp e e e B A 1 



F vx {t b )={ 



x 



(P + 2) -, 



Jy, for 2 < p, 



0.84 x 10 a x (1.64X 10 2 ) p ~ 2 > 
(S-fcJfl+^f- 1 ^ 



x f p e e AT 



-2 k 



2-k 



for 1 < p < 2, 



(17) 

where we have kept all model and observable parame- 
ters in the expressions. The flux density calculation has 
included a correction factor ps 0.2 due to both the ra- 
dial internal structure of the blast wave and the prolate 
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TABLE 2 

The sub-sample that is consistent with our coasting-in-wind model and their model parameter constraints. 



GRB 


ti 

do 2 s") 


h 

(TO 3 s") 


do 4 s~i 


F v X (*(>) 
(TO -7 Jvl 


z 


V fc e fc £j /max 


21 15 
e,-l t S,-2/ 


A 

7 17 

Cin- 2 e 22 e 22 ) 


14 1 
(-1Q52 11 11 \ 
V ±u t e, — l t B,-2 ei °/ 


050319 


61.1 


11.2 


8.5 


9.3 


3.24 


109.4 


49.5 


3.1 


7.6 


050416A 


2.5 


1.7 


26.2 


23.8 


0.65 


3.5 


45.1 


7.7 


5.1 


050713B 


7.9 


10.8 


47.9 


30.3 




8.3 


49.9 


8.9 


30.1 


050726 


4.2 


1.2 


1.7 


104.0 




3.1 


80.7 


2.6 


6.7 


050801 


0.7 


0.3 


4.6 


69.9 




0.7 


74.0 


2.0 


0.8 


050822 


64.1 


67.0 


52.3 


2.8 




168.5 


31.3 


9.0 


29.7 


051016B 


47.8 


66.4 


15.0 


1.1 


0.94 


188.5 


26.1 


8.3 


20.2 


060109 


7.4 


4.9 


1.8 


24.3 




9.2 


56.7 


3.8 


9.7 


060204B 


40.6 


5.5 


9.9 


26.8 




47.8 


54.9 


4.8 


12.3 


060210 


39.0 


24.2 


86.2 


18.3 


3.91 


47.8 


48.7 


7.1 


30.2 


060306 


2.5 


4.7 


12.4 


22.8 




3.2 


52.9 


4.6 


8.6 


060428A 


2.3 


11.0 


27.1 


27.7 




2.5 


50.9 


7.6 


28.5 


060507 


30.0 


7.0 


8.6 


8.7 




57.6 


45.5 


3.8 


5.8 


060526 


10.9 


10.0 


32.3 


6.1 


3.21 


23.9 


42.4 


3.9 


1.6 


060604 


35.2 


11.4 


40.4 


6.2 


2.68 


75.5 


40.4 


1.9 


6.2 


060707 


53.2 


22.2 


81.4 


1.9 


3.43 


186.0 


32.4 


4.5 


3.8 


060708 


38.1 


6.7 


43.9 


8.4 




74.6 


40.9 


5.6 


5.3 


060714 


3.2 


3.7 


33.2 


29.8 


2.71 


3.7 


55.3 


4.8 


6.8 


060719 


2.8 


9.6 


18.2 


9.2 




5.1 


43.4 


5.1 


8.9 


060729 


4.2 


73.0 


222.1 


23.6 


0.54 


4.1 


34.8 


42.1 


447.9 


060804 


1.8 


0.9 


12.2 


133.7 




1.2 


75.2 


4.2 


5.8 


060805A 


2.3 


1.3 


7.6 


5.9 




6.1 


45.0 


2.1 


0.6 


060813 


0.9 


1.8 


7.4 


225.1 




0.1 


83.0 


5.1 


22.0 


060814 


5.7 


17.4 


39.9 


20.7 




6.8 


46.6 


9.0 


37.6 


060912 


4.2 


1.1 


8.7 


31.6 


0.94 


5.3 


53.4 


1.6 


3.4 


070129 


13.2 


20.1 


54.6 


6.9 




25.5 


37.8 


7.9 


16.5 



Note. — The subscript "max" means the upper limit and "min" the lower limit. In deriving those constraints, we assumed a common 
value set p = 2.4 and k = 2. Observed red shifts are taken from Liang et al. (2007); for those with unknown red shift, z = 2 is assumed. 
A Hq = 71, £7a = 0.73, Q.m= 0.27 universe is assumed. 



equal-arrival-time surface effects (Granot, Piran & Sari 
f999). 

In principle, one can calculate these predictions and 
work out constraints using all the available parameters 
[p, k, F Vx (tb), etc.]. However, the appearance of p and k 
in the exponents would then obscure the dependence on 
basic parameters such as e e and e^, and make compar- 
isons very difficult. We choose instead to adopt the most 
common set of values k — 2 and p = 2.4 (corresponding 
to (3x,i — 1-2). Table [5] lists the individual constrain- 
ing results for the consistent sub-sample using individual 
values of t\, tb, F Vx (tb) and z. 

As an example for demonstration, in the following we 
use all typical observable values, i.e., p = 2.4, k = 2, 
h = 500 s, t b = 10 4 s, h = 10 5 s, v x = 1 keV, z = 2, and 
F Vx {tb) — 1 /^Jy, to carry out the constraints. It turns 
out the results using common values are consistent with 
those listed in Table [5] using individual values. Adopting 
these values, the items to be constrained become 



E iso = 0.5 x 10 44 A^t b erg. 



(18) 



v m {t x = 500s) = 5.7 x 10 15 ele 1 ^ 2 Al ,2 Tl Hz < v Xl (19) 



v c {t 2 = f0 5 s) = f .1 x 1 9 e B 3/2 A r 3/2 Tl Hz < is x , (20) 



F ux (t b = I0 4 s) = 5.9 x 10" 



7/5 1/10 ,ll/10 r 22/5 t 



Eq. flU) gives 

A* = 7.0 x I0 3 e ; 14 /ii eB i/" 
Plugging the above into Eq. (fT^|) , we have 



r -4 
1 o ■ 



,15/11 5/11 



< 0.9. 



(21) 



(22) 



(23) 



Various data modeling work and numerical experiments 
gave various values of e e and tB- A fairly conservative 
range is e e ~ 0.01 — 0.5, and es ~ 0.0001 — 0.1, respec- 
tively. Thus, the above constraint of e e and es can be 
easily met. Plugging Eq. (|22[) into Eq. (|2"0"|) gives an 
upper limit on Tq: 



r < 46 e, 



-21/88 15/88 
,-1 e B,-2- 



(24) 



When plugged back into Eq. (1221) , the above implies 

A, > 4.6 x 10- 2 e;J_/ 22 e B % 22 . (25) 
Finally, plugging Eq. (|22|) into Eq. dl8|) . we have 

E iso = 1.0 x 10 53 e~i 4 / n es er §- ( 26 ) 
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Fig. 5. — The total energy release in prompt 7-rays vs. the total 
energy of the ejecta inferred from X-ray data for our consistent 
sub-sample, both in isotropic equivalent values. Filled circles - 
bursts with measured z. Open circles — bursts with unknown z, for 



which z = 2 is assumed. The dotted line iS E*y 



For properties of wind from a typical Wolf-Rayet star, 
the lower limit on A* in Eq. (|25[) is easily met. Therefore, 
the major model parameter constraints from the above 
are T < 50 and E iso ~ 10 53 erg. 

We plot the individual E- lso vs. -E 7i i so — the isotropic 
energy release in prompt 7-rays — for the consistent 
sub-sample in Figure [5l An almost linear correlation 
between the two can be seen, which is consistent with 
and likely derives from the observed correlation between 
the plateau X-ray fluence and the prompt 7-ray flucncc 
(Liang ct al. 2007). It also shows that -E 7i i SO is compa- 
rable to or smaller than Ei SO , which alleviates the trou- 
blesome issue one faces in refreshed shock models of ex- 
tremely high 7-ray radiative efficiency ~ 90% (Zhang et 
al. 2007) . 

For all bursts in the sample, the post-plateau light 
curve does not show any further steepening break that 
can be identified as the jet break up to t%. This fact sets 
a lower limit on the ejecta beaming angle 6j > l/r^) = 

10 5 

21/88 -15/88 



r r 7 1 (t2Ab) (3 " fe)/[2(4 " fc)1 . For k = 2, tt « 10 4 s, t a 
s and using Eq. ([2"3]). we obtain 9j > 0.04 e e _ x c B 
radians — not a strong constraint. 

5. SUMMARY AND DISCUSSION 

By analyzing the closure relations for a large sample of 
X-ray afterglows with plateaus, we find that the plateau 
feature and post-plateau decay can be explained by the 
coasting-in-wind model in 55% of the sample. This sim- 
ple model is also able to explain the contemporaneously 
observed optical afterglow emission when optical data are 
available; however, for a very few bursts for which the end 
of the plateau is chromatic, additional assumptions (such 
as evolution of the microphysical parameters e e and es, 
or a two-component outflow) are needed to explain the 
lack of a break in the optical light curve. For the remain- 
ing 45% of the sample, the coasting-in-wind model can 
still reproduce the plateau, but the post-plateau decay in 
these bursts is too rapid to be explained in this model; a 



refreshed shock remains the most capable interpretation. 

Given the constraints derived from the consistent sub- 
sample, the coasting-in-wind model poses two physical 
challenges which must be addressed. 

First, is it reasonable for the wind medium to extend to 

the large distances r = 4r 2 ,c(W 2 ) 1/2 /( 1 + z) 10 18 cm, 
implied by this model? In fact it is: for a wind speed 
10 a v Wj s cms -1 , the wind ram pressure at this distance is 
10 -83 A*i;^ 8 (10 18 cm/r) 2 dynecm -2 , which exceeds the 
hydrostatic pressure (P ~ G£ 2 ) for any column den- 
sity S < 0.3AI /2 v wS (10 18 cm/r)gcur 2 . For the low- 
est acceptable values of A* we might therefore expect 
the wind to have terminated, especially within the high 
pressures and column densities of a starburst environ- 
ment. However, a fiducial wind (A* ~ 1) is sufficient to 
compete with starburst pressures (note that Dai & Wu 
2003 found A* ~ 0.01 for a burst which lacked an X- 
ray plateau). In fact, it is quite possible that at these 
radii the circum-burst medium is a merged wind from an 
entire star cluster (Chevalier & Clegg 1985), as massive 
stars are rarely found alone. The cluster mass function 
(McKee & Williams 1997) implies that the progenitor 
star is equally likely to belong to a massive cluster as a 
small one, and star clusters — especially massive ones 
— are very effective at clearing their immediate molec- 
ular environments with winds, radiation pressure, and 
photoevaporation before any stellar cores collapse (e.g., 
Krumholz & Matzner 2009). 

One must also consider the expansion of the wind- 
termination shock as the bubble expands (see Ramirez- 
Ruiz ct al. 2005); this is > 10 18 cm at 10 4 years if ci- 
ther uh < 10 2 j4*u* , g cm -3 in a uniform medium, or if 
the bubble expands into a previous stellar wind which is 
< 10 2 - 7 A*v^8 times denser (Koo & McKee 1992). How- 
ever, other cluster stars help to clear the ambient medium 
and alleviate this constraint. The duration of phase 2 is 
therefore not a strong constraint on the model, except in 
cases where there is independent evidence that the wind 
is weak and the ambient pressure and density are very 
high. 

The second and more puzzling challenge involves the 
arrangement of ejecta from the central engine. Our upper 
limit To < 50 is well below the lower limit of 10 2 ~ 3 de- 
rived from the prompt 7-rays using the pair opacity con- 
straint (Lithwick & Sari 2001) for a few Fermi Large Area 
Telescope (LAT) bursts, e.g., T > 600 - 900 for 080916C 
(Abdo et al. 2009a), T > 10 3 for short-hard burst 090510 
(Ackermann et al. 2010), T > 10 3 for 090902B (Abdo et 
al. 2009b), and T s» 200 - 700 for 090926A (Ackermann 
et al. 2011). However, the above constraints assume 
that both the MeV and GeV emission are produced in 
the same region; considering two-zone models for these 
two spectral components reduces the lower limit on T by 
a factor of ~ 5 (Zhao, Li & Bai 2011; Zou, Fan & Piran 
2011; Hascoet et al. 2011). 

The discrepancy on the T constraint between the above 
individual results and ours may be resolved by two 
means: (1) a possible selection bias. Three of the above 
four Fermi-LAT bursts do not have early (t < 5 x 10 4 
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s) X-ray or optical data to tell whether a plateau fea- 
ture is present. The last one, GRB 090510, has X-ray 
and optical coverage from t ~ 100 s and both X-ray 
and optical light curves show a shallow-to-steep transi- 
tion at t ~ 10 3 s which, however, is found to be con- 
sistent with a jet break (Kumar & Barniol Duran 2009). 
These four bursts are all among those most energetic ones 
(E- 



10 54-55 ergS; eyen the s hort-hard GRB 090510 
7, iso — 10 53 ergs). It is quite plausible that these 



^7, iso 

has E, ■ - 1 " 53 



most energetic bursts have much higher outflow LFs, so 
that their plateau features are too short (tb oc -EisoFq 4 ; 
see Eq. [TJ]) to be caught or identified. Future accu- 
mulation of Fermi-LAT bursts with early X-ray / opti- 
cal coverage could either support or disprove this bias 
effect. (2) The above constraints on T for each burst 
are derived from individual pulse(s) during the prompt 
phase, while our work is for the afterglow phase. It is 
very natural to have an ejecta bulk LF in the afterglow 
phase that is much lower than that of an 'individual emit- 
ter' during the prompt phase. This can be understood 
in the frames of two major GRB prompt phase models. 
In the conventional internal shock model (e.g., Paczyn- 
ski & Xu 1994; Rees & Meszaros 1994), discrete shells 
with large LF variation among them collide with each 
other, with each collision corresponding to an individual 
prompt pulse. Approximately after the prompt phase, all 
shells merge together and external shocks develop. It is 
possible that the prompt emission is produced only from 
the high LF portion of the outflow, while the outflow to- 
tal energy is dominated by the lower LF portion so that 
external shock possesses a lower LF. In another category 
of models that involve 'individual emitters' due to, e.g., 
magnetic turbulence inside a relativistic bulk flow (Lyu- 
tikov & Blandford 2003; Narayan & Kumar 2009), the 
prompt-emission-derived LF limit or value is actually for 
the bulk flow LF times the 'individual emitters' LF. Af- 
ter the 'individual emitters' energy is released in form of 
7-rays, the external shock LF is just that of the initial 
bulk flow. 

Alternatively the bulk LF may have been reduced by 
baryon pollution during time between the prompt and 
afterglow phases. Such contamination cannot come from 
the wind, for this would imply that the fast ejecta de- 
celerate early and skip the observed plateau. Possible 
stellar sources include: (1) the blowout of an inflated co- 
coon within the stellar envelope, which may have mixed 
with an unknown amount of stellar envelope material, 
and (2) the cap of stellar envelope which is unable to es- 
cape the advancing jet head, and is cast forward by the 
process of runaway shock acceleration (Matzner 2003). 
Of these, it is unclear how the former would pose an ob- 
stacle to the high-r jet material after it has escaped the 
star. 

The latter source, a trapped portion of the outer en- 
velope, is worthy of closer inspection. To evaluate this, 
we examine the approximate conditions for trapping of 
matter in front of the jet head (eqs. [26] and [27] of 
Matzner 2003, subject to his eq. [60). When trapping- 
occurs while the head's motion is nonrelativistic, the re- 

3 Note the typo (inverted expression) in the middle expression 
of equation (6) in Matzner (2003). 



suit is simple: the jet traps material once it is within 
~ 2.56jR* of the stellar surface, where i?* is the star's 
radius; in a polytropic region with index n, the trapped 
mass scales as 6* 3+ ™ (a very strong function of 6j). If 
instead the head is relativistic, the dependence weakens 
because a narrower, faster jet traps matter from deeper 
in the envelope; however it also depends on the intensity 
as well as the opening angle of the jet. For an example, 
consider the outer envelope profile for the SN 1998bw 
progenitor by Woosley, Eastman, & Schmidt (1999) and 
discussed by Matzner (2003, §5.1.1). In this compact, 
luminous, helium-stripped Wolf-Rayet star the isotropic 
rest energy of trapped material is 



3° 



1.0 



2.5 



10 52 ergs" 



0.63 



x 10 49 erg. 



(27) 

Although this material is swept into an accelerating for- 
ward shock (Matzner 2003), a negligible fraction attains 
LFs higher than Tj. Moreover, although it is small, 
the rest energy is well above the critical value r~ 2 £\ so 
required to decelerate the ejecta. Because it is about 
10~ 2 (10 52 erg/i?i sc ,) of the jet energy (for the fiducial case 
cited), the final LF after the interaction would be ~ 10 2 , 
but strongly dependent on 9j] our inferred Tj ~ 40 is 
entirely plausible as a final value. We refer the reader 
to the discussion by Thompson (2006), who considers a 
similar scenario. 

Whatever the origin of the lower- Tj matter, it is clear 
that it inherits effectively all of the GRB kinetic energy 
from the outflow powering the prompt phase. First, Fig- 
ure [5] shows that the afterglow energy budget is compat- 
ible with what is expected from the prompt phase, given 
reasonable radiative efficiencies. Second, if any signif- 
icant fraction of the energy had proceeded beyond the 
coasting shell at higher LF, it would have decelerated 
early and caused a noticeable departure from the plateau 
phase at early times. Indeed, since in the deceleration 
phase, when max(i/ m , i/ c ) < i>x (from Eqs. 1141 and [T7"|) . 



E iso o,F Vx {t b f /{p+2) tf P ' 2)/{2+P \ 



(28) 

we infer that the persistence of the plateau from times 
t\ to th rules out any early injection of fast ejecta with 
more than a fraction (t x I t h ) (?> ^) I ^v) ^ q.03 of the 
total energy. 

For 45% of the our sample, the post-plateau decay 
cannot be explained by phase 2 of the coasting-in-wind 
model. If we invoke the refreshed shock model for these 
afterglows, the inferred parameter suggests that the to- 
tal outflow kinetic energy is dominated by either the 
late ejecta (this issue might be alleviated for a lepton- 
dominated late ejecta) or by the slower, massive ejecta. 
This poses a great theoretical challenge to how the cen- 
tral engine works, and in fact is also the major moti- 
vation for us to investigate the coasting-in-wind model 
which does not face this serious constraint. The stan- 
dard rate history of the accretion due to the fallback 
of supernova-shocked stellar material gives a fairly steep 
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slope: M oc t~ 5 / 3 (Michel 1988; Chevalier 1989; Mac- 
Fadyen, Woosley & Heger 2001), which is unable to pro- 
vide conditions for a refreshed shock. It was recently 
proposed that an extended central engine activity could 
be realized due to the continuing accretion of the entire 
progenitor stellar envelope (Kumar, Narayan & Johnson 
2008a, b) or the slow self-adjustment of a transient debris 
torus formed at the stellar radius (~ 10 11 cm; Cannizzo 
& Gehrels 2009). For the second variety of the refreshed 
shock scenario where a late engine activity is not needed, 
the slower, massive ejecta that catches up with the de- 
celerating blastwave might actually be the broken-out 
cocoon that was produced when the GRB jet penetrat- 
ing through the star (Ramirez-Ruiz, Celotti & Rees 2002; 
Matzner 2003). 



Note added in proof: an independent investigation of 
the same physical model, with application to two bursts, 
has been posted in a preprint by Lei et al. (2011). 
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Enrico Ramirez-Ruiz and Eleonora Troja for useful com- 
ments about the manuscript. The research of R.-F.S. 
and C.D.M. is supported by an NSERC Discovery grant. 
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